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ABSTRACT 

Estimating the state of the power system in Real time is a ehallenging task, mainly beeause, solving the non-linear 
large power system equations by iterative solution eomsumes huge eomputational time and memory. Online Energy 
Management System (EMS) these state estimators need to Estimate the state of the system at very short intervals. 
Conventional methods like Newton-Raphson is not suitable for sueh real time applieation. 

To improve the performance of SE, eertain techniques like AI, Real time, Hirarehial and Dynamie methods are 
listed in reference [1][4][5][11][12][13]. The most eommon "Two-Level" SE technique is applied by spliting the large 
power system in to sub-networks and eaeh sub-networks are managed by the loeal eontrol station in eoordination with the 
Centre station - [2] [3] [7] [8] [9] [14]. Even though it reduees the eomputational time, still it is not suitable for real time 
applieations. In the "Two-Level SE" mainly same NR method is used, however, before applying the NR technique, the 
eleetrie network is physieally divided, whieh results to mathematieal approximation. 

This paper presents an approaeh to solve the problem in hand without physieally dividing the network and/or any 
approximation on the existing NR method. Single large mathematieal problem ean be divided into smaller independent 
parallel task by ehanging the multiplieation method of large sparee jacobian matrix. This leads to a whole new approaeh 
resulting in redueing in eomputing time as well as dynamie memory requirement. 

KEYWORDS: SE-State Estimation, WLS-Weight Least Square, RTU-Remote Terminal Unit, NR-Newton Rapson, 
AI-Artificial intelligenee, Node Area-A Node along with its Connected Node is Referred as Node Area 

INTRODUCTION 

The power system is highly distributed in nature, Instead of dividing the network in to sub networks, it is better to 
divide the large single problem into small independent sub tasks so that multiple proeessor ean be used for eomputation. 
Major ehallenge here to obtaing the independent sub-task. The NR method mathematieal solution is proeedure oriented 
iterative technique whieh involvs large number of non-linear simultaneous equations. The ideal solution would have been 
the vertical division of the probem into many sub-independent parallel task. By modifying the existing Newton-Raphson 
state estimation eomputational proeedure requires the elear understanding of the proeedural steps involed in the existing 
NR method. The brief insight of NR method is given below. 

NEWTON RAPSON (WLS) METHOD FOR STATE ESTIMATION -BRIEF REVIEW 

The Newton Rapson final equation is A Ax = b (taking weighting matrix W whieh is a diagonal matrix Wj; is the 
standard deviation of meter i.) By applying the tylor series to the nonlinear equations of power syatem following equations 
are derived. Here, the error vector ineludes nonlinear vector function f(x). In order to estimate x, an initial value x 0 is 
assumed and after every iteration x 0 will be up-dated till the convergence. 
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(J 0 T WJ)Ax = J 0 T WAz 



Let A= (J 0 T W J) & b = J 0 T W Az 



A Ax = b 



(1) 



This is a set of linear equations, if higher order terms of the taylor expansion of f(x) were really negligible, the 
solution yield the eorreet 'x'. The jacobian J is itself a function of x. The state variable vector x ean be obtained by solving 
the equation A*Ax = b iteratively. The vector x should therefore be ehanged aeeordingly after every iteration till the 
convergence is obtained. 

x n+1 = x n + Ax' (l st iteration eount) until convergence is reaehed. Weighted least squares (WLS) minimization 
technique is used. 

No of state variables = (2*n-l) 

Total no of measurements = m 

These measurements may inelude one or all quantities sueh as [ P ; ,Qj ,py ,qy, Vj ,5J 

• P; ,Qi = Real & Imaginary part of injected power respectively 

• py ,qy = Real & imaginary part of line follows respectively 

• V; & 5 ; = Voltage magnitude & phase angles. 

• Dimension of Jacobian matrix = m * (2n - 1) 

a measured ealeulated 

Az = z — z 

'W' is the diagonal weigh matrix of the order of (m*m) 

X= \5 7 ,d~,....d 7,v 7 v, v ; state variables. 

L 1 2 n — 1 1, 2, n,i 

n= number of network nodes, 1 ,2, . . .n. 

Elements of jacobian are derived from injected power and line flow equations. 
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PROPOSED ALGORITHUM 
Proposed New Algorithum 

The following modification of forming the jacobian results in independent parallel task and thereby multiple 
proeessors ean be used to earry out the task parallelly. Jacobian being the large sparse matrix, the suggested new method 
ean fully avoid storing and retrieving of large saprse matrix. the let 'm' be the total number of measurements taken for the 
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given network. Eaeh row of jacobian is formed based on the type of measurement undereonsideration, if there are 'm' 
measurements, the jacobian of the network will eontain 'm' rows and (2n-l) eoloms. 

Instead of forming the eomplete jacobian of the netork and then multiplying to fmd 'A' matrix, (see eql), It ean 
be eomputed by eaeh row of jacobean [J L *Wj */j = Aj ]and [ /j *Wj * Zj — b x ], this results in first terms of eaeh 
elements of the resultant matrix "A" and "b" respectively, whieh ean be visualized as shown below. 

{A 1+ A 2 ...+Ai...+ A m } Ax={b, ..+b r ...+ b m } (3) 



H H 

A= {A, + A 2 +A r .+ A m } 
b= {b,.+b r + b m } 



(4) 

(5) 
(6) 



Henee, it is not neeessary to form eomplete two dimentonal Jacobian and then earry-out the multiplieation task. 
New jacobian is a one dimention arry of the order of [l*(2n-l)] is replaeed by [m*(2n-l)], thereby no speeial sparse 
matries storage is required. 

Node Wise Grouping of Measurement: Jacobian Relation 

The measurements ean be grouped node wise, and it ean be represented as shown below. 



Node — l Measurements 
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A =UA P + A Q + ZA pij +Z A™ +A v + A S 
NAr \\ r r r r r r 



= J T *w * j 

J NAr "r NAr 



h =)\ b P +b Q +i: b pij + Y.b qij +b v +b 5 

NAr \\ r r r r r r 



(10) 



J NAr *W r * Z NAr (11) 

A = \ A NAl + A NA2 + - A NAr + •• + A NAn ) = X ( A NAj ) 

7=1 

* = + + - + *>AAr + - + ^NAn ) = I (*>AA/ ) 

7=1 

£(A^-)Ax= E(*/YA/) (12) 
7=1 7=1 



(13) 



(14) 



E(A ; ) = I(A M7 ) = A 

7=1 7=1 
m n 

7=1 7=1 
Henee AAx = b 

( -subseript 'NAi'refers toithnodearea 
p 

- It is the sub set of matrix 'A' for the iryeeted Real power measurement taken at 'r'th node/bus. Similarly, 
(here Q represents reactive power) and like wise for other measurements. 

Henee, from the above derivation, it is evedient that up to the formation of 'A NAi ' (node level), eomputations ean 
be earried-out independently without any approximation, whieh gives the seope for parallel proeessing. 

System Design 
Assumptions 

• All measurements are taken at the same instant 

• At eaeh node total of measurement ean be: - injected power measurements of node + (node voltage measurement 
+ node angle measurement +the line flow measurements ) of that node to other eonneeted nodes. 
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Flow Chart 
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Figure 1 

COMPUTING TIME ANALYSIS OF SE 
General Technique 

The eomputational time of state estimation using NR method is direetly linked with the number of nodes in the 
network. Generally, to solve A*Ax =b, if the eomputational time required is (say) "t," for the network having n^ number of 
nodes, then the eomputational time required for the network having n 2 number of nodes is approximately equals to 
(n 2 /ni) 2 *t!. Or in other words the SE eomputational time inereases "square time" the ratio of inerease in the number of 
nodes of the network. It's beeause as the number of node inereases the size of the Jacobian and "A" matrix. For example 

The size of matrix "A" for the network having 10 nodes is nearly equal to (20 x20) =400. 

The size of matrix "A" for the network having 30 nodes is nearly equal to (60 x60) =3600. 

(3600/400) = 9 = (30/10) 2 . 

Optimum Computing Technique 

The fact is, these matriees are highly sparse in nature. These none zero elements of jacobian/"A" has a direet 
relation with the network ineident matrix. Using this information one ean design the eomputation focusing only on none 
zero elements to avoid operations whieh results in zero. This ean be done by using the network ineident matrix as index to 
eompute arithmetie operation of none-zero elements. This makes the eomputational time required is direetly proportional 



52 



H. Nagaraja Udupa & H. Ravishankar Kamath 



to the number of measurements and does not depend on the size of the matriees. SE eomputational time reduees drastieally 
by using optimum eomputing technique. 

For example, a 10 bus system having number of lines eonneeted (say) = 14, then the size of Ybus = 10x10 = 100, 
whereas none zero elements in the Ybus = 10+14x2 = 38; 

Henee, none zero elements in the matrix "A" = 38 x4=152; whereas size of "A" approximately 
equal to (20x20) =400. Say "ti" is SE eomputational time for the network having 10 nodes using general method, then by 
using optimum technique, SE eomputational time for the same network is approximately equal to (152/400)x tj , let this be 
equal to t 2 . It's obvious that t 2 « t,. 

Now if the network size is inereased to n 2 (say for example n 2 =30, number of eonneeted lines=41, none zero Ybus 
elements=30+41x2=112 henee number of none zero elements of matrices= 112x4=448), then the eomputing time is 
approximately proportional to (448/152) = whieh is nearly equal '3', henee the eomputing time is nearly equal to= (n 2 /ni) x 
t 2 . Beeause, whatever the size of matriees, eomputational time depends upon the number of none zero elements whieh 
inereases almost linearly with node. 

Total SE eomputational time ean be related with the following equation 

{SE) t = (Axi = A~ l *b\ + f^J T */) + (s\ 

[SE) t = totalSE eomputational time 

\Axi = A -1 *b\ = Time taken to solve Axi 

fj T *Jj = Time taken to obtain A & b 

(s ) t = updating time 

EXAMPLE & RESULTS 

Note: All quantities are in pu. 

Table 1: Input Line Data for 13 Bus System (ISE) 









1-2 


.00147967 


.00286760 


2-3 


.00043843 


.00124174 


3-4 


.00027711 


.00078428 


4-5 


.00059760 


.00166769 


4-8 


.00159967 


.00310017 


5-6 


.00034314 


.00097190 


5-9 


.00034314 


.00097190 


L 6 " 7 


.00032364 


.00091669 


7-10 


.00032364 


.00091669 


8-9 


.00029438 


.00083380 


9-10 


.00053157 


.00150562 


9-12 


.00037793 


.00107050 


10-11 


.00058777 


.00166488 


11-13 


.00032364 


.00091669 


12-13 


.00036843 


.00104355 



Table 2: Injected Power & Voltage Measurements 



Ems 


Real 


lmg. 


Bus 
No. 


Volt 


1 


28.5335 


13.0837 


1. 


1.05300 
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Table 2: Contd. 



2 


-4.73 


-1.55 


3. 


0.95792 


3 


-1.27 


-0.41 


10. 


0.92304 


4 


-0.35 


-0.11 


11. 


0.91939 


5 


-4.38 


-1.44 


12. 


0.92178 


6 


-2.11 


-0.69 


13. 


0.92036 


7 


-0,42 


-0.13 


0 


0 


8 


-4.73 


-1.55 


0 


0 


9 


-1.27 


-0.41 


0 


0 


10 


-0.35 


-0.11 


0 


0 


11 


-4.38 


-1.44 


0 


0 


12 


-2.11 


-0.69 


0 


0 


13 


-0.42 


-0.13 


0 


0 



Table 3: Line Flow Measurements 



1-2 


28.5327 


13.0837 


2-1 


-27.2179 


-10.5355 


2-3 


22.4879 


8.9855 


3-2 


-22.2197 


-8.2258 


3-4 


20.9000 


7.7700 


[ 4-5 


12.9700 


5.1000 


4-8 


7.6200 


2.1282 


5-9 


4.3979 


1.9040 


L 8 " 4 


-7.4182 


-1.9161 


L 8 " 9 


2.7900 


0.3800 


9-10 


1..8043 


0.5194 


11-13 


-1.4600 


-0.4582 



Table 4: SE Result of L3 Bus Test System 
*Slack Bus 



No 




1 


1.053269 


0.101912 


2 


0.979368 


0.041320 


3 


0.958231 


0.015757 


4* 


0.945912 


0.000000 


5 


0.928889 


—0.021096 


6 


0,925904 


—0.025075 


7 


0.924513 


-0.026839 


8 


0.926468 


—0.022792 


9 


0.925278 


—0.025320 


10 


0.923398 


—0.028199 


11 


0.919737 


—0.033261 


12 


0.922133 


—0.029788 


13 


0.920707 


—0.031856 



Convergence =.001 for voltage, 



= .00001 for angle 



Table 5: Computational Time for 13 Bus Test System 



Conventional method 


1.08 ms 


0.88 ms 


New method (single proeessore) 


0.52 ms 


0.32 ms 


New method (paralle proeessores)* 


0.225 ms 


0.025 ms 



*Number of parallel proeessors are assumed to be =13 
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Table 6: Relative Time 





Total SE 


Time Ratio- 




Time 


Jacobian 




Ratio 


Produet 


Conventional /new single ~ 


2 


2.75 


Conventional /new parallel ~ 


5 


35 


New single/New parallel ~ 


2.36 


13 



The above table shows aetual results obtained whieh is same as that of the Integrated SE results obtained by the 
conventional method. It shows that both approaeh yields the same results. Thereby this new eoneept has been proved with 
this example. 

CONCLUSIONS 

Node wise grouping of data leads to new approaeh for applying the NR solution. The number of state variables to 
be eomputed at eaeh bus-level is very small as eompared to the whole network, and till the last step eomputations ean be 
earriedout independently & separately thereby the fast eomputation is possible. If sufficient numbers of measurements are 
ensured at eaeh node, then state estimantion ean be earried out at eaeh bus. At Central eomputing station only the final 
estimated state variables are to be eonsidered as sudeo-measurements. The other possibilities of this approaeh is under 
progress. 
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